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In this talk results from a new QCD analysis in leading (LO) and next-to- 
leading (NLO) order are presented. New parametrizations of the polarized 
quark and gluon densities are derived together with parametrizations of 
their fully correlated la error bands. Furthermore the value of a s (M§) 
is determined. Finally a number of low moments of the polarized parton 
densities are compared with results from lattice simulations. All details of 
the analysis are given in Ref. [1]. 



1. Formalism 

The twist-2 contributions to the structure function g\[x, Q 2 ) can be repre- 
sented in terms of a Mellin convolution of the polarized singlet density AS, 
the polarized gluon density AG, the polarized non-singlet density Aq , j vs , 
and the corresponding polarized Wilson coefficient functions AC- 4 by 

gi(x,Q 2 ) = 




where ej denotes the charge of the jth. quark flavor, Nf the number of 
flavors, and jXt is the factorization scale. In addition the above quantities 
depend on the renormalization scale \i r of the strong coupling constant 
a s (^) = ^)/(167r 2 ). 

The change of the parton distributions w.r.t. the factorization scale 
fJL 2 = Q 2 is described by the evolution equations which contain the polarized 
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splitting functions APjj. Both the polarized Wilson coefficient [2] and the 
polarized splitting functions [3] are known in the MS scheme up to NLO. 

The evolution equations are solved in Mellin-A space. A Mellin- 
transformation turns the MELLIN convolution into an ordinary product. 
The solutions for the evolved parton distributions are structured such that 
the input and the evolution part factorize [4], which is of key importance for 
the error calculation. An inverse MELLiN-transformation to x-space is then 
performed numerically by a contour integral in the complex plane around 
all singularities. 

2. Parametrization and Error Calculation 

The shape chosen for the parameterization of the polarized parton distribu- 
tions at the input scale of Q 2 = 4.0 GeV 2 is: 

xA qi (x,Ql) = 7]iAiX a *(l -x) h (l +-fiX + piX^). (2) 

The normalization constant A} is defined such that 77, is the first moment 
of Agj(x, Qo)- The densities to be fitted are Au v , Ad v , Aq, and AG. 

Assuming SU(3) flavor symmetry the first moments of Au v and Ad v 
can be fixed by the SU(3) parameters F and D measured in neutron and 
hyperon /3-decays to r] Uv = 0.926 and rjd v = —0.341. The sea-quark distribu- 
tion Aq was assumed to be described according to SU(3) flavor symmetry. 
Given the present accuracy of the data we set a number of parameters to 
zero, namely p Uv = p dv = 0, 7 5 = pg = 0, and j G = p G = 0. Further- 
more we adopted the following two parameter relations: a G = (iq + c, with 
0.5 < c < 1.0 and (bq/bc)(pol) = (bq/bG){unpol). These relations were 
essential to respect positivity for Aq and AG and to achieve the expected 
similar low-x behaviour for both parton densities. No positivity constraint 
was assumed for Au v and Ad v . In addition Aqcd was determined. 

The gradients of the polarized parton densities w.r.t. the fitted param- 
eters needed for the Gaussian error propagation are calculated analytically 
at the input scale Qq in Mellin-A^ space. Their values at Q 2 are then 
given by evolution. For more details see Ref. [1]. 

To treat all data sets on the same footing only the statistical errors were 
used. To be able to calculate the fully correlated la error bands via Gaussian 
error propagation only fits ending with a positive-definite covariance matrix 
were accepted. We allowed for a relative normalization shift between the 
different data sets within the normalization uncertainties quoted by the 
experiments. Thereby the main systematic uncertainties of the data were 
taken into account. The relative normalization shifts being fitted once and 
then fixed enter as an additional term (penalty) for each data set in the 
X 2 -expression. 
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3. Results 

The results reported here are based on 435 data points of asymmetry data, 
i.e. gx/F\ or Ai, above Q 2 = 1.0 GeV 2 , the world statistics published so far 
(see Ref. [1] for a list of references). The fits are performed on g\ which is 
evaluated from the asymmetry data using parametrizations for Fi [5] and R 
[6]. The data do not constrain the four parameters j Uv , 7^, bq, and be well 
enough. Their values have been fixed as obtained in the first minimization. 
The NLO polarized parton densities at the input scale are presented in Fig. 
1 (ISET = 3, see Ref. [1]). 
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Fig. 1. Polarized parton distributions at the input scale Q\ = 4.0 GeV 2 (solid line) 
compared to results obtained by GRSV (dashed-dotted line) [7] , AAC (dashed line) 
[8]. The shaded areas represent the fully correlated \a error bands calculated by 
Gaussian error propagation [1]. The dark dotted lines correspond to the positivity 
bounds given by the unpolarized densities. 

The current data constrain Au v at best, followed by Ad v , Aq, and AG. 
Mainly the lack of data at low x leads to the broader error bands for the 
latter two densities. The comparison with results from other QCD analyses 
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shows the variation of the outcoming parton distributions when using differ- 
ent parametrizations at the input scale. The measured structure function g\ 
is well described both as function of x and Q 2 . The derived parton densities 
and their error bands have been evolved to Q 2 values up to 10,000 GeV 2 . 
One finds that Au v stays positive, Ad v remains negative although less con- 
straint than Au v , and AG stays positive even within the error band for the 
range Q 2 > A GeV 2 . All three densities evolve towards smaller values of 
x. The sea-quark density Aq is negative for Q 2 < 10 4 GeV 2 and remains 
negative within errors for x < 5 • 10~ 2 up to Q 2 = W 4 GeV 2 , but changes 
sign for larger values of x. 

In determining a s the QCD-parameter Aqcd was fitted. The impact of 
the variation of both the renormalization and the factorization scale on the 
value of a s was investigated. The following result was obtained (for ISET 
= 3, see Ref. [1], as well for more details): 

/»^ o-,-,o +0.004 . . +0.004 ,, , +0.008 , , 
a s (M z ) = 0.113 _ Q 004 (stat) _ Q Q04 (fac) _ 005 ( ren ) ■ 



This value is compatible within la with the world average of 0.118 ± 0.002 
[9], although our central value is somewhat lower. It is also compatible with 
results from other QCD analyses of polarized and unpolarized data. 

In recent lattice simulations [10] low moments for the polarized par- 
ton densities Au v , Ad v , and Au — Ad were determined. In Table 1 these 
moments are compared with the ones extracted from our NLO polarized 
densities. 







QCD results 


lattice results 
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NLO moments 


QCDSF 


LHPC/ 






at Q 2 = 4 GeV 2 




SESAM 


Au v 


-1 


0.926 ± 0.071 


0.889(29) 


0.860(69) 







0.163 + 0.014 


0.198(8) 


0.242(22) 




1 


0.055 ± 0.006 


0.041(9) 


0.116(42) 


Ad v 


-1 


-0.341 + 0.123 


-0.236(27) 


-0.171(43) 







-0.047 + 0.021 


-0.048(3) 


-0.029(13) 




1 


-0.015 + 0.009 


-0.028(2) 


0.001(25) 


Au-Ad 


-1 


1.267 + 0.142 


1.14(3) 


1.031(81) 







0.210 + 0.025 


0.246(9) 


0.271(25) 




1 


0.070 + 0.011 


0.069(9) 


0.115(49) 



Table 1: Moments of the NLO parton densities (for ISET =3, see [1]) at Q 2 = 
4 GeV 2 and from recent lattice simulations at the scale ji 2 = 1/a 2 ~ 4 GeV 2 . 



proceed printed on February 1, 2008 



■5 



The first moment is denoted with n = — 1. For the n = 0, 1 values of 
the QCDSF Collaboration no continuum extrapolation was performed. The 
values compare within errors. Unlike in the unpolarized case where a strong 
m^-dependence is expected in the lattice extrapolation the comparison here 
suggests a flat behaviour, which also has been found in [11] very recently. 

4. Conclusions 

A QCD analysis in LO and NLO of the current world-data on polarized 
structure functions was performed. New parametrizations of the polarized 
parton densities including their fully correlated la error bands were de- 
rived. These parameterizations are available as a fast FORTRAN-routine which 
makes their application possible in Monte Carlo simulations. The value de- 
termined for a s (Mz) is compatible within la with the world average. Com- 
paring the lowest moments with values from lattice simulations the errors 
improved during recent years and the values became closer. 
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